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Abstract 

The production of black holes in extra-dimensional brane-world theories can lead to 
detectable signals via the Hawking evaporation of the black hole to brane-localised Stan- 
dard Model modes. We calculate, as a function of partial wave number and number of 
toroidally compactified extra dimensions, the leading correction to the energy spectrum of 
such Hawking radiation (the greybody factors) for decay into spin 1/2 fermions and spin 
1 gauge fields localised on the Standard Model brane. We derive the associated improved 
differential emission rates for both types of fields. We provide both simple expressions 
for the leading behavior of the greybody factors in the low-energy limit ujth <^ 1, and 
numerical evaluation of our full analytic expressions for the emission rates, valid for 
uoth ~ 1- The full analytic expressions demonstrate that both the greybody factors and 
emission rates are enhanced as the number of extra dimensions increases. 



1 Introduction 



One of the most spectacular consequences of brane- world theories [1, 2, 3] with low 
fundamental gravitational scale, > TeV (for some works on the structure of these 
theories, on their cosmological implications and experimental signatures, see [4]), is the 
modification of the properties of small black holes [5] that, in principle, allows them to be 
copiously produced at the LHC and in other high-energy processes [6, 7, 8, 9, 10, 11, 12]. 
In a collider setting, the dominant signal for such black hole production arises from the 
Hawking evaporation products of the black hole, in particular the characteristic spec- 
trum of such radiation. In a previous paper [13], we considered the leading semiclassical 
corrections to the spectrum of Hawking radiation for a higher-dimensional black hole 
evaporating to scalar modes. We considered both a scalar field localised to the brane as 
well as the situation in which it could propagate in the full (4 -|-n) -dimensional bulk, and 
computed the corresponding greybody factors as functions of energy o;, angular momen- 
tum j, and number of (flat) extra dimensions n. The scalar case is both a useful testing 
place of the methods used in such a computation, and potentially relevant for Higgs, or 
other scalar boson production from black holes formed in collisions at the LHC or other 
high energy colliders. However, from a practical point of view the most important cases 
to consider are those of the production of SM fermions or SM gauge bosons to which we 
turn in this paper. 

In both cases the fields that need to be considered are brane locahsed. For fermions 
this is because of the chiral nature of the SM fermions, since a fermion propagating in the 
bulk is necessarily non-chiral from the 4d perspective.^ It is possible to consider SM gauge 
fields propagating in the higher-dimensional bulk, or a subspace thereof [15]. However, 
the already existing collider constraints bound the size of such gauge extra dimensions 
to -Rgauge ~ 3 TeV, while, in the best possible case where only the graviton propagates 
in the extra dimensions, the limits on the new fundamental higher-dimensional Planck 
scale are of order > 1.2 TeV [16]. On the other hand, the validity of the semiclassical 
description of black holes requires that their horizon radius be large compared to the 
higher dimensional Planck length th ^ 1/M*, which given the just cited bounds, as well 
as the hmits on the maximum mass of black hole that can be produced at the LHC, 
implies that rn must be somewhat large compared to i^gauge; too. In this case the gauge 
fields can also be considered as effectively brane localised, the 'brane' in question having 
some non-trivial substructure. To leading order this substructure does not change the 
greybody factors we compute for the spin 1 case. 

Note that the properties of a black hole, in particular its Schwarzschild radius, th, 
and production cross section, differ from the 4d result provided that rn is smaller than 
the size of at least some of the extra dimensions rn < R- If the brane tension does not 
strongly perturb the (4 -|- n)-dimensional black hole solution, and the extra dimensional 

^The RHD neutrino can propagate in the bulk, leading to an alternate, non-seesaw explanation of 
the lightness of the observed neutrino states [14]. Thus in principle it is possible to consider the Hawking 
decay of a black hole to these bulk neutrino states. However, from the perspective of the brane observer 
this leads only to a missing energy signal which is not distinguishable from the many other sources of 
missing energy in brane- world theories. Therefore we do not consider this possibility in the present work. 
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space is approximately flat on the scale r^, then the black hole is well-described as a 
angularly symmetric (4-|-n)-dimensional black hole centered on the brane, but extending 
out into the n extra dimensions, leading to the expression [5, 17], 

for the Schwarzschild radius. Here, we follow these commonly made assumptions, which 
have the consequence of limiting our discussion to the Hawking decay of black holes 
in the context of Arkani-Hamed-Dimopoulos-Dvali like theories [1]. In particular, black 
holes in theories of the Randall- Sundrum type [3] require a more involved treatment (for 
a recent paper considering the decay of RS black holes see [18]). 

After such black holes are produced, they decay by the emission of Hawking radiation, 
and it is expected that they will decay mostly to particles on our brane [19]. The spectrum 
of radiation, with a temperature Tbh = (1 + ^)/(47rrii-) {in G = = c = h = 1 units), 
is given by the Hawking formula 



ir-j^, ^^'"^"^^ exp {u:/Tsh) T 1 WY^' ' ^ ^ 



where j labels the total angular momentum quantum number, h labels any other quantum 
numbers of the emitted particle, as well as the particle type, and in the phase-space 
integral = a; for a massless particle. Here a{uj), known as the 'greybody factor', is 
an energy-dependent function arising from the backscattering of part of the outgoing 
radiation back into the black hole by the non-trivial metric in the region exterior to the 
black hole horizon (for early works on the emission spectra of 4d black holes, see [21]). 

In this paper we calculate the greybody factors and emission rates for spin 1/2 and 
spin 1 particles localised to the SM brane in the presence of a non-rotating black hole 
of the type described above. It is expected that such a non-rotating black hole is a good 
approximation to the exact solution after the initial 'balding' and 'spin-down' phases [7]. 
As discussed for the scalar case in our first paper [13], the greybody factor is computed 
using its equality with the absorption cross section for the appropriate type of particle 
incident on the background metric that describes the brane black hole. The equality of 
greybody factors to absorption cross sections implies that the greybody factors do not 
invalidate the thermal nature of the black hole. 

Finally, it is important to recall that the semiclassical calculation of Hawking emission 
is only reliable when the energy of the emitted particle is small compared to the black 
hole mass uj <C M, since only in this case is it correct to neglect the back reaction of the 
metric during the emission process. This in turn requires that the Hawking temperature 
obeys the relation Tbh ^ M ^ which is equivalent to demanding that the black hole mass 
M 3> M*. Inevitably this condition breaks down during the final stages of the decay 
process, but for those black holes of initial mass larger than M* most of the evaporation 
process is within the semi-classical regime. 

In Section 2 of this paper, we present the general equation for the radial part of the 
wave function of fields with spin 1/2 and 1, as well as a discussion of the metric used in 
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calculating the greybody factors. We proceed to solve this equation at the asymptotic 
regimes of the near-horizon and far-field zones, for arbitrary number of extra dimensions 
and angular momentum number, and we match them in an intermediate zone in order 
to construct a complete, smooth solution over the whole radial regime. 

In section 3, the greybody factors for fermionic and gauge fields are computed. Our 
full analytic results are contained in Eqs. (35-38), (40), (46), and (49). Although analyt- 
ically complicated, these full results are necessary for the accurate determination of the 
Hawking emission spectrum from branc black holes (see section 4). We expect they will 
be of use to those searching for black holes at current and future colliders, or in ultra- 
high-energy-cosmic-ray collisions. However, for ease of comparison with previous studies 
of greybody factors we derive simplified expressions for the leading corrections to the 
spectrum of the Hawking radiation valid for urn <^ 1- These expressions are presented 
in Eq. (47) for spin 1/2 and Eq. (55) for spin 1, for arbitrary j and n. Their values for 
n = 2,4,6 extra dimensions, and for the first three, and thus most important, angular 
momentum modes (with j = i, |, | for spin 1/2, and j =1, 2, 3 for spin 1) are tabulated 
in Table 1 and Table 2, respectively. 

In section 4, we proceed to calculate the full greybody factors and emission rates 
for spin 1/2 and spin 1 fields by a numerical evaluation of our full analytic results of 
sections 2 and 3 (valid beyond the simplified limit uth ^ 1). A similar analysis for spin 
fields, based on the results derived in [13], is also included in this section. The results 
are presented in figures 1,2 and 3. Our conclusions are summarized in section 5. 



2 General equation for spin 1/2 and 1 fields and its 
solution 



The purpose of this paper is to calculate the greybody factors for the emission of fermion 
and gauge fields in the background of the metric on their 4-dimensional brane induced 
from a (4 + n)-dimensional Schwarzschild-like black hole centered on the brane. The 
starting point of our analysis will be the general equation for the propagation of fields 
with spin s — 1/2 and 1 in an arbitrary 4-dimensional spherically-symmetric background. 
Using the Newman- Penrose formalism [22] (see [23] for a compact review), this general 
field equation may be written as 



Following the analysis of Cvetic and Larscn [23], it is possible to derive from this general 
equation a set of 2s differential equations for the components of the wave function of the 
spin s field. Specifically, assuming a general, spherically-symmetric background of the 
form 



AB1...B2S-1 



0. 



(3) 




(4) 
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and using a factorized ansatz for each component of the wave function of the field of the 
form 

= A-'/^ (2S)'=/2- Sk-si^) , (5) 

the set of coupled differential equations separates into two sets of decoupled equations for 
the radial and angular parts of the wave function. In the above expression, k — 0, ...,2s 
labels the helicity of the components, and S^si^) — e^^'^ d^^^(6) , where j and m are 
the total angular momentum number and its projection on a fixed axis, respectively, and 
^im stand for the rotation matrices. Focusing on the radial part of the wave function, the 
form of the corresponding equation turns out to depend on the helicity of the component. 
For example, for the upper component ijjs {k = 2s), it takes the form [23]: 



dr 



dr 



— + 2isujdr^-A + 



In the above, A = j {j + 1) — s {s — 1). As we will shortly see, solving the above equation 
for the upper component of the field only is sufficient to lead to the determination of the 
absorption cross section, and thus the greybody factor. 

In the presence of additional space-like dimensions, the higher-dimensional spheri- 
cally symmetric generalizations of the 4-dimensional Schwarzschild solution have a line 
element given by 

ds^ = -h{r) df + h{r)-^ dr^ + dQl+„, (7) 

where 

= (8) 
The angular part of the above higher-dimensional metric tensor is 

dQl^^ = del^^ + sin' On+i (del + sin' ( - + sin' 62 {dOl + sin' dif^) ...)], (9) 



with < if < 2tt and < < tt, for i = 1, ...,n + 1. It is easy to check that, at 
distances r ^ R, where R the size of the extra compact dimensions, the above metric 
tensor reduces to the usual 4-dimensional Schwarzschild solution. 

Under the assumption that the Standard Model fields, both fermions and gauge 
bosons, are restricted to live on a 4-dimensional brane, they propagate in a gravitational 
background which is given by the induced 4d metric on the brane following from the 
higher-dimensional black-hole background (7). This follows if we set On — t^/2, forn > 2. 
Specifically it takes the form 

Comparing Eqs. (4) and (10), wc may easily see that the general metric functions ap- 
pearing in Eq. (4) are now given by: E = r' and A — hr^. Substituting in the general. 
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radial differential equation (6), we obtain the following simplified equation: 



2\l-s 



dr 



,2^2 



+ \ + 2isur — (n + 1) (1 - /i) - 



-A - (2s - 1) (s - 1) (n + 1) (1 - i?, = . (11) 
In the above, we have used the relation 

f=^-^{l-h), (12) 
dr r 

that follows from the definition Eq. (8). In this paper, we are going to study the cases 
of fermion (s = 1/2) and gauge-boson fields (s = 1). For those fields, the last term 
appearing inside the curly brackets in Eq. (11) vanishes independently of the number of 
projected extra dimensions and, therefore, it will be ignored in the following analysis. 

In order to compute the grcybody factor for the Hawking radiation in the gravita- 
tional background described in Eq. (10), through the emission of fermion and gauge-boson 
fields, we need to solve Eq. (11) for the radial part of the emitted field in the region out- 
side the horizon of the black hole and all the way to infinity. Due to the complexity of 
the differential equation, the exact solution appears impossible to derive. As in [24] [13], 
(see also [25] for related works) we are going to follow an approximation method that 
is suitable for low energies urn ^ 1- The method involves solving Eq. (11) first in the 
vicinity of the black hole, then at infinite distance from it and, finally, matching the two 
asymptotic solutions in an intermediate regime. This method will help us derive the form 
of the absorption coefficient in the equivalent scattering problem of an incoming wave 
propagating in the background (10), in terms of which the absorption cross section, and 
thus the greybody factor of the emitted radiation is defined. 

As we will shortly see, different components of the emitted field carry different parts 
of the particle wave function. For an emitted field with spin s 7^ 0, it is only the upper, ipg, 
and lower, ijj-s, components that are the radiative ones. The upper component will turn 
out to correspond mainly to the incoming wave while the lower component corresponds 
mainly to the reflected wave. Moreover, the equation for the lower component is the 
charge conjugate of the one for the upper component, therefore solving for any one of 
them provides the solution for the other one. Here, we will concentrate on the incoming 
wave and compute the absorption coefficient by comparing the incoming fiux at the 
asymptotic regimes of the horizon and infinity. The matching of the two solutions in the 
intermediate regime, however, is still necessary. We impose the following normalizations 
for the incoming mode of a field, with spin s, 

at the horizon of the black hole, where Pbh — ^/Tbh, and 

^(00) ^ ^(00) ^2uj^Y''^ ^-i'^r ^^4^ 

at infinity. In the forthcoming sections, and for the solution of the general, radial equation 
(11), we closely follow the analysis presented in Ref. [13] for the case of the emission of 
scalar fields. 
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2.1 Solving the general equation in the near- horizon zone 



In the vicinity of the black hole, Eq. (11) can be solved by making a change of variable, 
r — > /i. Using Eqs. (8) and (12), we may write the general-spin field equation in the form 



\ /-. , \ + 2s) , 1 dR 
^ ^ ^ ^ (n+l) \ dh 



IS UJTh 



{n+iyh{l- h) (n + 1)2(1 -/i) {n+l)h 



R^O, (15) 



where, for simplicity of the notation, the subscript s has been dropped from R. By 
using the redefinition R{h) = /i"(l — h)'^F{h), the above equation takes the form of a 
hypergeometric equation 



r/^F dF 
h{l-h)-— + [c-{l + a + b)h] — -abF^O, 
dn'^ dh 



with 



a — a -\- (5 -\- 



s + n (1 — s) 



h — a-\-(5, c—l — s-\-2a. 



(16) 



(17) 



(rz + 1) ' 

The power coefficients a and /3, in turn, are found by solving second-order algebraic 
equations leading to the results 



a^ — s-\- 



lUjrH 
n + 1 ' 



"n + 1 ' 



and 



1 - 2s ± ^(1 + 2j)2 - Aujh-jj - SisujrH 



2{n+l) 

respectively. The general solution of the hypergeometric equation (16) is 



(18) 
(19) 



RNH{h) = {1-hf F{a, b, c; h)+A+ /i"" (l-h)^ F{a-c+l, b-c+1, 2-c; h) , (20) 

where A±. arc arbitrary constants. Expanding the above solution in the limit r ^ th 
(equivalently /i ^ 0) we obtain the asymptotic behaviour 



for a — a+, and 

Rnh — A_ exp(-ia;r^+^y) + A+ cxp (^zcjr^+^y) 
for a — a^. In the above, we have used the 'tortoise' coordinate y defined as 

hi //(/•) 
^"rr(r^ + l)' 



(21) 



(22) 



(23) 



The choice a = a+ leads to a solution with an outgoing wave of zero amplitude at the 
horizon, and an incoming wave with infinite amplitude. This is an irregular solution which 
must be discarded. On the other hand, the choice oi a — leads to regular incoming 
and outgoing waves with amplitude unity at the horizon. The boundary condition that 
only incoming modes are to be found in the region outside the horizon of a black hole 
forces us also to set = 0. After these choices are made, the solution near the horizon 
has exactly the normalization of the fields described in Eq. (13), since q;_ = —iuj Pbh/^t^i 
thus determining the asymptotic normahzation constant to be A^^"* — A_. Turning to the 
(5 coefficient, the criterion for the convergence of the hypergeometric function F(a, 6, c; h) 

Re (c - a - 6) = - ( ^(1 + 2j)2 - Auj'^rjj - Siscurn ) > (24) 

clearly demands that we choose /3 = 

We finally need to shift the hypergeometric function towards large values of r. This 
may be done, by using a standard linear transformation formula [26], in the following 
way: 



RNHih) = A_ /i" 



r(l-. + 2a)r(-2/3 + i-^) 



^1 - hf ^7 "^'^ , F(a, b,a + b-c+l;l-h) 

' ^ r(a-/3 + l-.)r(a-/3+i-f) ^ ' ' ' ^ 



r(l - s + 2a)T(2(3 



l-2s 



+ (1 - /,)-Wi-2.)/(n+i) ^ V ^^^^ y s+ lTiBk ^i'^ -a,c-b,c-a-b + l;l-h) 



r(a + /?)r(a + /?+ 



(25) 



We can now "stretch" the above expression towards the intermediate regime, by expand- 
ing in the limit r — > oo, or /i — > 1, and take 

RMk) . A_(^)-^^-^ r(i-. + 2.)r( 2, + ^) 

r(«-/3 + l-.)r(a-/3+i^) 



,+1 r(l-s + 2a)r(2/3 



Note that, in order to simphfy the procedure of the "matching" , the low energy hmit 
has been taken in the expression of the f3 coefficient in the power of r. No expansion has 
been made, up to this point, in the arguments of the Gamma functions. 



2.2 Solving the general equation in the far-field zone 

We now need to find the far-field zone solution before being able to match the two 
asymptotics in the intermediate regime. Going back to the general equation (11), and in 
the limit r ^ oo or /i — > 1, we obtain 

(PR 2{l -s)dR /2 2isa; A\^_q 
dr'^ r dr \ r r^J 
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In order to solve the above equation, we set: R = e^*"^*^' r-'^* R{r). By also making a change 
of variable z = 2ia;r, the above equation adopts the form of a confluent hypergeometric 
equation 

(PR /, X dR ~ , , 

+ (6 - ^) — - ai? = , (28) 



dz^ dz 
with a — j — s + 1 and b — 2j + 2, and general solution 

R{z) = B+ M{a, b, z) + S_ C/(a, 6, z) , (29) 

where M and U are the Kummer functions, and B± are arbitrary coefficients. We may 
therefore write the complete solution for the radial function at infinity as 

Rppir) = e-'^'^'r^^' [s+M(j-s + l, 2j + 2, 2iujr)+B_U{j - s + l, 2j + 2, 2iu;r)] . (30) 

We first need to expand the above expression for large values of r. This expansion 
will help us make sure that the emitted fields have the far-field normalization given in 
Eq. (14) and will enable us to compute the incoming fiux at infinity. In the limit r — > oo, 
we find 



icLir „2s— 1 



s+l 



B_ + 



B+ e'^^^-'+^^T{2j + 2) 



T{j + s + l) 



^iur B+T{2j + 2) 

^ T{j-s + l) {2iujy+'+^r ^ 

(31) 

The first part of the above expression gives the incoming wave at infinity, and thus, by 
comparison with Eq. (14), defines the coefficient A^°°\ The second part gives the outgoing 
wave. This outgoing part is suppressed as it should be. We remind the reader that the 
above solution corresponds only to the upper component of the emitted field (cither to 
ijji/2 or to ijji) which in turn corresponds to the incoming wave. The charge conjugate 
of the above solution gives the lower component (the other radiative component) of the 
emitted field, which will have a dominant outgoing wave and a suppressed incoming one. 

We may now stretch the far-field solution (30) towards small values of r. In the limit 
a;r ^ 1, we therefore obtain 

«-W-a.^'^' + ;^f(-^|^^|^- (32) 

Matching the two solutions (26) and (32) in the intermediate regime, we obtain the 
relations 

B = A_ r(i-. + 2a)r(-2/? + i^) 

ri+^r(a-/3+l-5)r(a-/3+i^)' ^ ' 

^ ^ A- ry^' {2tLoy^^' r(l - . + 2a)T{2P - igf ) T{j -s + l) 

r(a + /5)r(a + /? + ^±^)r(2j + i) 

With the determination of the above arbitrary coefficients, we have now completed the 
determination of the complete solution for both fermions and gauge-boson fields prop- 
agating in the background of a projected (4 -|- n)-dimensional Schwarzschild black hole 
on a 4-dimensional brane. We now proceed to the calculation of the greybody factor in 
each case for an arbitrary number of extra compact space-like dimensions. 
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3 Greybody factors for the emission of fermionic and 
gauge fields 



As we will shortly see, the expression of the greybody factors for Hawking radiation 
for fermions and gauge fields, involves the quantity determined by the two 

normalization coefficients appearing in the asymptotic solutions (13) and (14). By using 
the relations (33)- (34), this takes the form 



AW 



A_ (2a; 



\3+s piTT (j-s+l)/2 



{'2urH) 



j+s p-i-K (j-s+l)/2 



r(l-s + 2a) 



r(j + s + i) 

where the coefficients C and D stand for 

22,+i gm(.-i/2) p(^2/3 - 1^) r(j - s + 1) 

^ ^ r(a + /3)r(a + /3 + ^±^)r(2j^T7' 

r(2j + 2)r(-2/j + ^) 

r(a-/? + i-s)r(a-/? + i^)r(j + s + i) 

The measure squared of the ratio A^^^ jAS^^ can then be written as 

J^{h) 2 



(35) 



D 



AM 



|r(l - s + 2q;)|2 CC*(u;r/^)4j+2 + (a;ri/)2J+i {CD* + C*L>) + L>i:>* 



(36) 



(37) 



(38) 



The derivation of a final explicit result for the greybody factors or emission rates from 
the above expression requires the evaluation of Gamma functions for complex arguments. 
This task is a tedious procedure with the result being highly complicated and unillumi- 
nating. Accurate evaluations of the greybody factors and associated emission rates do 
require use of the above expressions, and so in section 4, we present numerical results 
for the final expression of the greybody factors, as a function of loth^ by using the exact 
value of the above ratio. However, /or comparison with previous studies, and as in the 
case of the analysis for the emission of scalar fields [13], a simplified expression may be 
derived in the low-energy limit uth ^ 1- Expanding the Gamma functions in this limit, 
wc may express the denominator of the fraction appearing in Eq. (38) in a power series 
in {ujru) keeping only the leading term. The result of the expansion depends strongly 
on the spin of the emitted particle, therefore, we will now discuss separately the results 
for each case. We emphasise that applications in which accurate values of the emissions 
rate are required should not use these simplified expressions valid only for urn <^ 1, but 
rather should use the full semiclassical results derived above (together with Eqs. (40), 
(46), and (49)) as in the evaluation presented in section 4. 



3.1 Spin 1/2 fields 

Following Ref. [23], we define the incoming flux of a fermionic field as the radial com- 
ponent of the conserved current, = y/2 a'^^ ip^ip^ , integrated over a two-dimensional 
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|^s=l/2| 



(40) 



sphere. This may be finally written as 

1 dN , ,o , ,9 , ^ 

5^^ = |ft/.r-|fl-V.r (39) 

in terms of the measure squared of the upper and lower components of the fermionic 
field. Taking the ratio of the value of this quantity at the horizon with respect to its 
value at infinity, we obtain the absorption probability for the scattering of an incoming 
wave travelling from spatial infinity towards the horizon of the projected black hole. It 
can be simply expressed as 

if we take into account the fact that the upper component contributes little to the infalling 
flux both at inflnity and at the horizon. 

Setting s = 1/2 and expanding the C and D coefficients, defined in Eqs. (36)-(37), 
in the limit luth -C 1, we obtain 

2-(4i+2)/(n+l) 

CC* = r(i + l)^ +0{u;rHr, (41) 
22j+i 

CD* + C*D = + 0{curHy, (42) 

TT 

04j 2(4j+2)/(n+l) 

DD* = T{j + lf + 0{ujrHf. (43) 

Taking into account the fact that the smallest physically allowed value of the total angular 
momentum j is the value of the spin, s = 1/2, we easily conclude that the dominant 
term in the denominator of Eq. (38), in the low-energy limit, is given by the first term 
in the expansion of DD*. Moreover, wc have 



cosh(2Q;7r) 1 2tt (oorH) .-^ , , ,s 

^ ^ ~ - + . \ ^' + Oiu^mf , (44) 



|r(l/2 + 2a)|2 TT TT ' (n+l)2 



where a is the imaginary part of a. Putting everything together, we obtain 
The absorption probability leads to the absorption cross section crabs(<^) via 

^abs ('^) = —2 I I (46) 

for any spin s and angular momentum partial wave number j. Thus, for the fermion case 
the final absorption cross section is, in leading order in {uirn), 

a/2). . 27rV|,(2i + l)(a;r^)^^-^2-(W/(n-n) 
^abs H = 2^^T(i + l)^ ■ ^ 
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This is the primary analytic result valid for corn <^ 1 for the case of spin 1/2 fermions. 

As it is obvious from the above formula, the absorption cross-section has a strong 
dependence on the total angular momentum number j of the emitted field as it affects 
both the power of the dominant (a;ri^)-term and the multiplicative coefficient in front. 
On the other hand, as in the case of the emission of brane-locahsed scalar fields studied 
in [13], no dependence of the power of (corH) on the number of extra dimensions emerges 
for a brane-localised fermionic field. Nevertheless, the multiplicative coefficient in front 
depends on the number of extra dimensions. 

In order to be able to draw concrete conclusions about the dependence of 
on j and n, the value of the absorption cross-section for n = 2, 4 and 6, and j = 1, 
I and I have been computed and displayed in Table 1. By looking at the entries of 
this Table for j — ^ and for every n, we see that the absorption cross section becomes 
proportional to the area of the horizon of the black hole, times a numerical coefficient, 
and is independent of the energy of the emitted particle. As the total angular momentum 
number j increases, the absorption cross section rapidly decreases as it gets suppressed 
by extra powers of urrn and its numerical coefficient also gets smaller. If instead we fix j 
and vary the number of extra dimensions, then we see that the absorption cross section 
gets enhanced as the number of the projected extra dimensions increases. The behaviour, 
therefore, of the leading term in the expression of the greybody factor, in the low-energy 
regime, for the emission of brane fermionic fields is similar to the one for the emission of 
brane scalar fields [13]. 



Table 1 : Absorption Cross Section for a brane — localised fermionic field 


n = 2 


j = 1/2 

J = 3/2 
J = 5/2 


^(1/2)^2^3^,2^ + ... 


n = 4 


J = 1/2 
j = 3/2 
j = 5/2 


aiif ^4 2V5^r|, + ... 
^iif ^ f 2^/^!, K^)^ + . . . 


n = 6 


J = 1/2 
j = 3/2 

J = 5/2 


^(V2)^4 23A^r|, + ... 

^iif ^|26/^7rr|,(u;rH)^ + ... 

ai;f)^^2Vv^rl,(c.rHr + ... 
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3.2 Spin 1 fields 

We now move to the case of gauge fields. In this case, no conserved current may be 

defined, however, the incoming flux can be expressed as the integral of the flow of energy, 

j^i j^f 

T^^^ — 2a'^^,a'^Q,ip^^ip , through a 2-dimensional sphere [23]. Then we obtain the 
expression, 

2^f = 2^0«'l"-|«-.r), (48) 

where the redshift of the local energy of photons has been taken into account. The 
absorption probability is given again by the ratio of the value of the above quantity at 
the horizon with respect to its value at infinity, and it is found to be 



\A 



s=l\ 



AW 



(49) 



Going back to Eq. (38) and setting s = 1, we expand the C and D coefficients in the 
limit luth <^ 1- Thus, we obtain 

cc* = + 0{u;rHr, (50) 



-i2^^+i sin(^)sin(^ 



(i+i) 



DD* = , .\? , J^2 + 0{ur^r. (52) 

The first term in the expansion of DD* gives again the dominant term in the denominator 
of Eq. (38) in the low-energy limit. Also, we have 



2a sinh(2Q;7r) 4 (cjr//) 



2 



OioorHT. (53) 



r(2Q;) r(-2Q;) tt (n + 1)2 

By using the above results, the absorption probability may be written as 

+ r(2gi)'r(2j + 2)- 

leading to the following absorption cross section, for the gauge field case, in leading order 
in ujrn, 

a« (.) - 4.r^ i2.rn?^ ( rfe) 1(^1) r(, + 2) \ ^ 
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This is our primary analytic result valid for urn <^ 1 for the case of gauge fields. 

As in the case of scalar and fermionic fields, the absorption cross section for brane 
localised gauge fields depends non-trivially on the angular momentum number j and the 
number of the extra dimensions n. The power of the dominant (a;r//)-term depends only 
on j but the numerical coefficient in front depends strongly on both parameters. The 
values of the above quantity for n —2, 4 and 6, and j —1,2 and 3 are given in Table 2. 



Table 2 : Absorption Cross Section for a brane — localised gauge field 


n = 2 


J = l 

i = 2 

J = 3 




n = 4 


J = l 

J = 2 
J = 3 


aW^A2V5r(^)2r(|)V|,(u;r^)^ + ... 

(1) r(3/5)2r(4/5)2 2 / , 
"abs — 39375 r(7/5)2 K^^H) + ■ ■ ■ 


n = 6 


i = i 

J = 2 

i = 3 





Unlike the case of scalar and fermionic field emission, the absorption cross section 
for the emission of gauge fields depends on the energy of the emitted particle even in the 
case of the lowest partial wave with j = 1. This follows from the different dependence 
of the dominant {uruYieim. on the angular momentum number. Any further increase in 
the value of j leads to the appearance of a suppression factor of powers of {wth) while 
the increase in the number of extra dimensions strongly enhances the absorption cross 
section. 

4 Evaluation of Results 

In this section, we proceed to calculate the corresponding emission rates for spin 1/2 and 
1 fields in the background of the projected 4-dimensional black hole given in Eq. (10). In 
this section we will not use the simplified expressions valid only for corn ^ 1, but rather 
the full analytic expressions Eqs. (36-38). 

For the relevant 4-dimensional process of the emission of brane-locahsed fields, Eq. (2) 
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a. 




2 



5 



4 



3 



1 
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Figure 1: The greybody factors a^hsi'-^) for the emission of (a) fermions, and (b) gauge 
brane-localized fields, as a function of corn, for n — 0,2,4 and 6. 

takes the form 



Thus for the evaluation of the emission rates, we need the expressions of the greybody 
factors summed over the angular momentum number j. A simple numerical analysis 
shows that, by summing over the first three partial waves, we obtain the dependence of 
the greybody factors on the parameter ujrH to high accuracy, with all higher partial waves 
adding an almost zero contribution to the final result. The numerical evaluation gives us 
the ability to use the exact value of the ratio (38) instead of the simplified leading-order 
corrections presented in Eqs. (47) and (55), which are valid only for ujrn -C 1. Such a 
numerical analysis will reveal how fast the low-energy approximation breaks down as the 
expansion parameter urn increases. Figures 1(a) and (b) depict the greybody factors for 
fermions and gauge fields, respectively, in units of 7rr|^, as functions of ujth- The first 
three partial waves have been summed in each case and the dependence for n = 0,2,4 
and 6 is shown. We observe that, for small enough values of tur//, we reproduce, as 
expected, the behaviour given by the leading-order corrections and according to which 
the greybody factors are enhanced as the number of extra dimensions projected on 
our brane increases. This agreement holds for both species of fields up to the value 
uoth ~ 0.4. As uoth increases further, it is the lower-dimensional models that seem to 
give the largest greybody factors. However, our calculation is limited to a WKB-like 
semiclassical approximation and only an exact numerical analysis of the original master 
equation, Eq. (11), can determine (Tabs(<^) in the intermediate and high-energy regime 
^fH ^ 1- Let us finally note that our analysis correctly reproduces the behaviour of the 
greybody factors in the case of n = [21, 27], with (T^j^{^''(c(j) adopting a non- vanishing 
asymptotic value at very low energies and c'"abs(<^) going to zero. The same behaviour is 
observed in the cases with n 7^ 0. 

Figures 2(a) and (b) display the behaviour of the differential emission rate per time 
unit and energy interval d'^E/dtdw versus the parameter corn, for fermions and gauge 
fields, as it follows by combining Eq. (56) with the exact expression of the greybody 
factor in each case and the definition of the temperature of the projected 4-dimensional 




(56) 
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0.1 0.2 0.3 0.4 0.5 0.6 q 0.2 0.4 0.6 0.8 

loth corn 



Figure 2: The differential emission rate (PE/dtdcu for (a) fermions and (b) gauge fields, 
as a function of corn, for n — 0,1,2, A and 6. 

black hole, Th = (n + l)/47rri^. The emission rates are substantially enhanced as the 
number of extra dimensions increases. The increase is larger for gauge fields which leads 
to the result that the emission rates for the two species become comparable, for large 
values of n. On the other hand, for low values of n, the emission rate for a fermion is 
considerably larger compared than that for a gauge field, with the difference reaching an 
order of magnitude in the limiting case of n = 0. 

It is useful to perform a similar numerical analysis for the case of scalar fields being 
emitted by the same 4-dimensional black hole, a case studied in Ref. [13]. By follow- 
ing a similar analytical approach, the leading-order correction in the expression of the 
greybody factor in that case was also derived. According to those results, the leading 
correction for the first, and thus dominant, partial wave was independent of the number 
of extra dimensions. Any difference in the value of the greybody factor was to come 
from the higher partial waves whose coefficient indeed increased with increasing n. As 
uoth increases, however, the next-to-leading-order corrections for each partial wave (for 
example, the terms being denoted by in Table 1 and 2, and similarly in Table II in 
[13]) can dominate the contribution coming from the higher partial waves. In the case of 




0.1 0.2 0.3 0.4 0.5 . 1 .2 .3 0.4 0.5 

torn ojrn 



Figure 3: (a) The greybody factor for the emission of hrane-localized scalar fields for 
n = 0, 1 and 6; (b) the corresponding differential emission rate for n = 0, 1, 2, 4 and 6. 
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fermions and gauge fields, we showed that this does not happen until the value loth ~ 
0.4. For scalar fields, it turns out that this value is much smaller. Figure 3(a) depicts 
the greybody factor for the emission of brane-localized scalar fields, by using the exact 
value of the ratio of the two coefficients and given in Eq. (46) of Ref. [13]. We 
immediately see that, for extremely small values of curH, the greybody factors, for various 
values of n, converge to the same asymptotic value, in agreement with our simplified 
expression. However, very soon, the next-to-leading-order corrections dominate and cause 
the greybody factor to decrease as n increases. The emission rate for scalar fields is 
nevertheless enhanced as the number of extra dimensions increases. This is caused by 
the decrease in the denominator of Eq. (56), due to the dependence of the temperature 
Tbh on n, that overcomes the decrease of the greybody factor. As in the case of fermion 
and gauge fields, the increase amounts to orders of magnitude compared to the limiting 
case of n = 0. 



5 Conclusions 

In an earlier paper [13], we studied the problem of the decay of a higher-dimensional 
Schwarzschild-like black hole through the emission of cither higher-dimensional bulk, 
or brane-localised, scalar modes. In this paper, we extended the computation of the 
greybody factors that appear in the Hawking decay process to the spin 1/2 and spin 1 
brane-localised cases of greatest phenomenological interest. Since the angular momentum 
of the black hole has been ignored, our calculations arc relevant to the post balding and 
spin-down phases of the life of a black hole produced at a high energy collider or by ultra 
high energy cosmic ray interactions. The semiclassical approximation methods employed 
allow us to calculate the greybody factors and the corresponding emission rates for a 
general number n of flat extra dimensions of the Arkani-Hamed, Dimopoulos, and Dvali 
type. (It is important to bear in mind that black holes in highly curved bulk spacctimcs, 
for example Randall-Sundrum 5d theories, require a somewhat modified treatment from 
the one presented here.) The primary analytic results of our calculation for the greybody 
factors, for arbitrary j and n, are presented in Eqs. (35-38), and (40), (46), and (49). 
Simplified expressions are provided in Eqs. (47) and (55) spin 1/2 and spin 1 respectively. 
In both cases there is an increase in the greybody factor as n is increased at fixed j. As 
the subsequent numerical analysis revealed, this behaviour survives up to intermediate 
values of the energy of the emitted particle and contributes to a considerably enhanced 
emission rate for both fermions and gauge bosons. A similar analysis for the emission 
of scalar fields also led to an enhancement of the emission rate as n increases despite 
the fact that the greybody factor, rapidly deviating from the behaviour dictated by the 
leading-order correction in the low-energy approximation, is suppressed. We expect that 
our results will be of use in the detailed study of the signatures of possible black hole 
production events. 

Note added. While writing this paper, the related work of Ref. [28] appeared which 
also considers the greybody factors for black holes in brane-world theories. 



16 



Acknowledgments. We wish to thank Roberto Emparan and Nemanja Kaloper for 
discussions. 



References 

[1] N. Arkani-Hamed, S. Dimopoulos and G. R. Dvah, Phys. Lett. B 429, 263 (1998) 
[hep-ph/9803315]; Phys. Rev. D 59, 086004 (1999) [hcp-ph/9807344]; 
I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, Phys. Lett. B 436, 
257 (1998) [hep-ph/9804398]. 

[2] K. Akama, Lect. Notes Phys. 176, 267 (1982) [hep-th/0001113]; 

V. A. Rubakov and M. E. Shaposhnikov, Phys. Lett. B 125, 139 (1983); Phys. Lett. 
B 125, 136 (1983); 

M. Visscr, Phys. Lett. B 159, 22 (1985) [hep-th/9910093]; 
1. Antoniadis, Phys. Lett. B 246, 377 (1990); 

1. Antoniadis, K. Benakh and M. Quiros, Phys. Lett. B 331, 313 (1994) [hep- 
ph/9403290]; 

J. D. Lykken, Phys. Rev. D 54, 3693 (1996) [hep-th/9603133]. 

[3] L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 3370 (1999) [hep-ph/9905221]; 
Phys. Rev. Lett. 83, 4690 (1999) [hcp-th/9906064]. 

[4] 1. Antoniadis, K. Benakh and M. Quiros, Phys. Lett. B 460, 176 (1999) [hep- 
ph/9905311]; 

N. Arkani-Hamed, S. Dimopoulos and J. March-Russell, Phys. Rev. D 63, 064020 
(2001) [hep-th/9809124]; 

R. Sundrum, Phys. Rev. D 59, 085009 (1999) [hep-ph/9805471]; 

G. F. Giudice, R. Rattazzi and J. D. Wells, Nucl. Phys. B 544, 3 (1999) [hep- 

ph/9811291]; 

E. A. MirabeUi, M. Perelstein and M. E. Peskin, Phys. Rev. Lett. 82, 2236 (1999) 
[hep-ph/9811337]; 

J. L. Hewett, Phys. Rev. Lett. 82, 4765 (1999) [hep-ph/9811356]; 

T. Han, J. D. Lykken and R. J. Zhang, Phys. Rev. D 59, 105006 (1999) [hep- 

ph/9811350]; 

N. Kaloper, J. March-Russell, G. D. Starkman and M. Trodden, Phys. Rev. Lett. 
85, 928 (2000) [hep-ph/0002001]; 

S. Cullcn and M. Perelstein, Phys. Rev. Lett. 83, 268 (1999) [hep-ph/9903422]; 
S. Hannestad and G. Raffeh, Phys. Rev. Lett. 87, 051301 (2001) [hep-ph/0103201]; 
L. J. HaU and D. R. Smith, Phys. Rev. D 60, 085008 (1999) [hep-ph/9904267]; 
N. Arkani-Hamed, S. Dimopoulos, N. Kaloper and J. March-Russell, Nucl. Phys. B 
567, 189 (2000) [hep-ph/9903224]; 

S. Cullen, M. Perelstein and M. E. Peskin, Phys. Rev. D 62, 055012 (2000) [hep- 

ph/0001166]; 

P. Kanti and K. A. Olive, Phys. Rev. D 60, 043502 (1999) [hep-ph/9903524]; Phys. 



17 



Lett. B 464, 192 (1999) [hep-ph/9906331]; 
A. Dedes and P. Kanti, hep-ph/0003051; 

R. Guedens, D. Clancy and A. R. Liddle, Phys. Rev. D 66, 043513 (2002) [astro- 
ph/0205149]; 

T. Matsuda, Phys. Rev. D 66, 107301 (2002) [hep-ph/0209214]; 
M. Bastero-Gil, V. Di Clemente and S. F. King, hep-ph/0211012. 

[5] P. C. Argyres, S. Dimopoulos and J. March-Russell, Phys. Lett. B 441, 96 (1998) 
[hep-th/9808138]. 

[6] T. Banks and W. Fischler, hep-th/9906038. 

[7] S. B. Giddings and S. Thomas, Phys. Rev. D 65, 056010 (2002) [hep-ph/0106219]. 

[8] S. Dimopoulos and G. Landsberg, Phys. Rev. Lett. 87, 161602 (2001) [hep- 
ph/0106295]. 

[9] S. Dimopoulos and R. Emparan, Phys. Lett. B 526, 393 (2002) [hep-ph/0108060]; 

G. F. Giudicc, R. Rattazzi and J. D. Wells, Nucl. Phys. B 630, 293 (2002) [hep- 

ph/0112161]; 

D. M. Eardley and S. B. Giddings, Phys. Rev. D 66, 044011 (2002) [gr-qc/0201034]. 

[10] M. B. Voloshin, Phys. Lett. B 518, 137 (2001) [hep-ph/0107119]; Phys. Lett. B 524, 
376 (2002) [hep-ph/0111099]. 

[11] S. Hossenfelder, S. Hofmann, M. Bleicher and H. Stocker, Phys. Rev. D 66, 101502 
(2002) [hep-ph/0109085]; Phys. Lett. 548, 73 (2002) [hep-ph/01 12186]; 

H. C. Kim, S. H. Moon and J. H. Yce, JHEP 0202, 046 (2002) [gr-qc/0108071]; 
K. Cheung, Phys. Rev. Lett. 88, 221602 (2002) [hep-ph/01 10 163]; Phys. Rev. D 66, 
036007 (2002) [hep-ph/0205033]; 

R. Casadio and B. Harms, Int. J. Mod. Phys. A 17, 4635 (2002) [hep-th/01 10255]; 
S. Hofmann, M. Bleicher, L. Gerland, S. Hossenfelder, S. Schwabe and H. Stocker, 

hep-ph/0111052; 

G. Pasztor and T. G. Rizzo, in Proc. of the APS/DPF/DPB Summer Study on the 
Future of Particle Physics (Snowmass 2001) ed. R. Davidson and C. Quigg, hep- 
ph/0112054; 

G. Landsberg, Phys. Rev. Lett. 88, 181801 (2002) [hep-ph/0112061]; 

E. J. Ahn, M. Cavagha and A. V. Olinto, Phys. Lett. B 551, 1 (2003) [hep- 
th/0201042]; 

U. Baur et al, in Proc. of the APS/DPF/DPB Summer Study on the Future of Par- 
ticle Physics (Snowmass 2001) ed. R. Davidson and C. Quigg, hep-ph/0201227; 
T. G. Rizzo, JHEP 0202, Oil (2002) [hep-ph/0201228]; 
S. N. Solodukhin, Phys. Lett. B 533, 153 (2002) [hep-ph/0201248]; 
V. Cardoso and J. P. Lemos, Phys. Lett. B 538, 1 (2002) [gr-qc/0202019]; gr- 
qc/0211094; hep-th/0212168. 

A. Chamblin and G. C. Nayak, Phys. Rev. D 66, 091901 (2002) [hep-ph/0206060]; 
V. Frolov and D. Stojkovic, Phys. Rev. D 66, 084002 (2002) [hep-th/0206046]; Phys. 



18 



Rev. Lett. 89, 151302 (2002) [hep-th/0208102]; 

S. C. Park and H. S. Song, hep-ph/0111069; 

Y. Uehara, hep-ph/0205068; hep-ph/0205199; 

A. V. Kotwal and C. Hays, hep-ph/0206055; 

T. Han, G. D. Kribs and B. McElrath, hep-ph/0207003. 

[12] A. Goyal, A. Gupta and N. Mahajan, Phys. Rev. D 63, 043003 (2001) [hep- 
ph/0005030]; 

J. L. Feng and A. D. Shapere, Phys. Rev. Lett. 88, 021303 (2002) [hep-ph/0109106]; 
L. Anchordoqui and H. Goldberg, Phys. Rev. D 65, 047502 (2002) [hep-ph/0109242]; 
R. Emparan, M. Masip and R. Rattazzi, Phys. Rev. D 65, 064023 (2002) [hep- 
ph/0109287]; 

L. A. Anchordoqui, J. L. Feng, H. Goldberg and A. D. Shapere, Phys. Rev. D 65, 
124027 (2002) [hep-ph/0112247]; 

Y. Uehara, Prog. Theor. Phys. 107, 621 (2002) [hep-ph/0110382]; 

J. Alvarez-Muniz, J. L. Feng, F. Halzen, T. Han and D. Hooper, Phys. Rev. D 65, 

124015 (2002) [hep-ph/0202081]; 

A. Ringwald and H. Tu, Phys. Lett. B 525, 135 (2002) [hep-ph/0111042]; 

M. Kowalski, A. Ringwald and H. Tu, Phys. Lett. B 529, 1 (2002) [hep-ph/0201139]; 

D. Kazanas and A. Nicolaidis, hep-ph/0109247; 

P. Jain, S. Kar, S. Panda and J. P. Ralston, hep-ph/0201232. 

[13] P. Kanti and J. March-Russell, Phys. Rev. D 66, 024023 (2002) [hep-ph/0203223]. 

[14] N. Arkani-Hamed, S. Dimopoulos, G. R. Dvali and J. March-Russell, Phys. Rev. D 
65, 024032 (2002) [hep-ph/9811448]; 

K. R. Dienes, E. Dudas and T. Gherghetta, Nucl. Phys. B 557, 25 (1999) [hep- 
ph/9811428]. 

[15] K. R. Dienes, E. Dudas and T. Gherghetta, Phys. Lett. B 436 (1998) 55 [hep- 
ph/9803466]; Nucl. Phys. B 537 (1999) 47 [hep-ph/9806292]. 

[16] See e.g., J. Hewett and J. March- Russell, Extra Dimensions Review, in Review of 
Particle Properties, Phys. Rev. D 66 (2002) 010001. 

[17] R. C. Myers and M. J. Perry, Annals Phys. 172, 304 (1986). 

[18] R. Emparan, J. Garcia-Belhdo and N. Kaloper, hep-th/0212132. 

[19] R. Emparan, G. T. Horowitz and R. C. Myers, Phys. Rev. Lett. 85, 499 (2000) 
[hep-th/0003118]. 

[20] S. W. Hawking, Commun. Math. Phys. 43, 199 (1975). 

[21] D. N. Page, Phys. Rev. D 13, 198 (1976); 
W. G. Unruh, Phys. Rev. D 14, 3251 (1976); 
N. Sanchez, Phys. Rev. D 18, 1030 (1978). 



19 



[22] S. Chandrasckhar, The Mathematical Theory of Black Holes (Oxford University 
Press, New York, 1983). 

[23] M. Cvetic and F. Larsen, Phys. Rev. D 57, 6297 (1998) [hep-th/9712118]. 

[24] S. R. Das and S. D. Mathur, Nucl. Phys. B 478, 561 (1996) [hep-th/9606185]; 
J. Maldacena and A. Strominger, Phys. Rev. D 55, 861 (1997) [hep-th/9609026]; 
S. R. Das, G. W. Gibbons and S. D. Mathur, Phys. Rev. Lett. 78, 417 (1997) 
[hep-th/9609052]. 

[25] S. S. Gubser, Phys. Rev. D 56, 7854 (1997) [hep-th/9706100]; 
R. Emparan, Nucl. Phys. B 516, 297 (1998)[hep-th/9706204]; 
V. Frolov and D. Stojkovic, gr-qc/0211055. 

[26] M. Abramowitz and I. Stegun, Handbook of Mathematical Functions (Academic, 
New York, 1966). 

[27] J. H. MacGibbon and B. R. Webber, Phys. Rev. D 41, 3052 (1990). 
[28] D. Ida, K. y. Oda and S. C. Park, hep-th/0212108. 



20 



